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ABSTRACT 

We present near infra-red (NIR) adaptive optics assisted spectroscopic observa- 
tions of the CO (A/i = 2) absorption bands towards the centre of the giant eUiptical 
galaxy NGC 1399. The observations were made with NAOS-CONICA (ESQ VLT) and 
have a FWHM resolution of O'.'IS (14pc). Kinematic analysis of the observations re- 
veals a decoupled core and strongly non-Gaussian line-of-sight velocity profiles (VPs) 
in the central 0.2 arcsec (19pc). NIR imaging also indicates an asymmetric elongation 
of the central isophotes in the same region. 

We use spherical orbit-superposition models to interpret the kinematics, using a 
set of orthogonal "eigenVPs" that allow us to fit models directly to spectra. The models 
require a central black hole of mass 1.2lQ'g x IO^Mq, with a strongly tangentially biased 
orbit distribution in the inner 40pc. 

Key words: instrumentation: adaptive optics, Galaxies: kinematics and dynamics, 
galaxies: individual: NGC 1399 



> 1 INTRODUCTION 



Super-massive black holes (SMBHs) are thought to be the 
only viable candidates for the massive dark object (MDO) 
observed at the centres of many nearby galaxies. Indeed, re- 
cent near infra-red (NIR) observations of the centre of the 
Milky Way have resolved individual stars orbiting in close 
proximity to the central MDO (which coincides with the ra- 
dio source, Sgr A*) and these data rule out all other plausible 
expla nations for a MDO, other than a SMBH iSchodel et alJ 

Significantly, a relationship between the mass of the 
SMBH and the bulge luminosity of the host galaxy was dis- 
covered ( Kormcndv & Richstonc 1995) and subsequently, a 
tighter correlation between the SMBH mass, M., and the 
velocity di spersion of the bulge, a (the M.-g r elation) was 
measured jFerrarese fc Merrilj|200ot lGebhardtl.2000) of the 
form 



log (M./Mq) = Q + /31og(cr/cro). 



(1) 



with cro=200km s~ MTremaine et alJ J200l hereafte r T02) 
find a = 8.13 and B = 4.02 and iFerrarese fc FordI ||200,'tI. 
hereafter FF05) find a = 8.22 and (3 = 4.86. Similar re- 
lations with low scatter have al so been found between M. 
and the infrared luminosity LTHj|Mgjggni_fc Hi|n j2n0^ and 
between M, and bulge mass tearing fc RiJ2004^ ■ 

It is believed that the mass accretion history of a 



SMBH is linked to the formation and evolution of its host 
(Hachnclt & Kauffmann 2000; dc Zccuw 2003) and so such 
a precise relation, connecting quantities on vastly different 
scales, provides an important constraint on models of galaxy 
assembly. It would be particularly remarkable if it holds 
true for galaxies of different morphological types, which 
most likely underwent ver y different fo r matio n and evolu- 
tion histories. For example, iFaber et al.l J^^^) suggest that 
power-law ellipticals and spiral bulges formed dissipatively 
whereas core-like ellipticals formed from mergers, yet both 
appear to follow the same relation. In practice, the rela- 
tion can also be used to measure the mass of a host galaxy 
black hole (BH) where a dynamical estimate is not possible 
jAller fc Ric.hstnnelEooilYrfc Tremain^lEnn^. 



1.1 Contention 

There has been co nsiderable debate over the values of the pa- 
rameters a and /3 ( Forr arese fc Mer rit 2000; Gcbhardt"20QlJ; 
iTremaine et~arji2002i:.Ferrarese fc Ford.2005.1 . The slope /? is 
crucial for comparison with theoretical models that attempt 
to explain the M,-cr relation, but currently, the sample of 
galaxies used by T02 (Fig. and FF05 is somewhat limited 
and biased. Of the 31 galaxies in the sample of T02, 18 are 
elliptical, 9 are lenticular and 4 are spiral (see Fig. Ql. Al- 
though there are roughly equal numbers of power-law and 
core ellipticals, there remain few high dispersion and low 
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dispersion galaxies where any deviations from the canonical 
slope will be most obvious. Such bias is not surprising con- 
sidering that, until recently, only one telescope could per- 
form dynamical mass estimates with the required spatial 
resolution: the Hubble Space Telescope (HST). 

In order to reliably estimate the mass of the black hole, 
it is necessary to resolve kinematics in the region of space 
where the black hole potential dominates over the potential 
of the stars, a point stressed by FF05. The radius of this 
sphere of influence (Sol) is of order 



where G is the gravitational constant and a* is the aver- 
age stellar velocity dispersion of the spheroidal component. 
Ground based observations have, in the past, been limited 
by atmospheric seeing to a resolution of around l"or more. 
There are very few galaxies close enough and with suffi- 
ciently large SMBHs to have a Sol larger than this, so nearly 
all of the sample of T02 and FF05 comes from HST data 
(diffraction limit of ~0'.'04 at 4000 A). However, HST has its 
limitations. 

Massive giant ellipticals which populate the high disper- 
sion region of the M.-cr plane are rare (so on average more 
distant) and have low su rface brightness c entres with flat 
core- like density profiles (|Faber et alJIlQQ^ . Although the 
total luminosity of such galaxies is high, they have too low a 
surface brightness to efficiently achieve high signal-to-noise 
(SNR) observations with HST's 2.4m primary. Conversely, 
high SNR nuclea r observations of g alaxies with steep power- 
law photometry l)Faber et al.lll99^ are feasible due to their 
large numbers (therefore proximity) and rapid increase in 
surface brightness towards the nucleus. Hence, less massive 
ellipticals with steep power-law surface brightness profiles 
are common in the sample of T02, but very massive ellipti- 
cals with core-like surface brightness profiles are rare. 

The low dispersion region of the M.-cr plane is popu- 
lated by low mass bulges in spirals that tend to be dusty and 
obscured at the centre. Extinction is reduced at longer wave- 
lengths (the extinction at 2.3fim is only 10% of that at V- 
band), but HST is unable to probe the nuclei of such galaxies 
as it does not have a suitable IR spectrograph. Therefore, 
less massive dusty spirals are also under-represented in the 
current sample. 

The outcome of these selection effects is an uneven 
distribution of galaxies along the a axis: 70% of the T02 
sample have 120km s~^< a < 250km s~^. It is impera- 
tive that galaxies with high and low velocity dispersions be 
investigated, both to verify the uniformity of the M.-cr re- 
lation over different morphological types and a larger dis- 
persion range, but also to better determine the nature of 
the relation. The slope of the M.-cr relation remains con- 
tentious (T02, FF05) because data points at the extremi- 
ties, crucial for defining the slope, are sparse. Furthermore, 
some attractive theories of galaxy evolution predict depar- 
tures from the current power-law relation which would only 
be d etectable with more data at high and low dispersions 
jHaehnclt & Kauffmann 2000; Zhao ct al. 2002). 

Modern adaptive optics (AO) facilities on large 8m class 
ground based telescopes (diffraction limit of 0'.'059 at 2.3fim) 
are the key to solving this problem. The spatial resolution 
achieved with adaptive optics is limited primarily by the 



diffraction limit of the telescope. In the case of NIR ob- 
servations with an 8m primary, this can match the spatial 
resolution of HST in the optical. Furthermore, 8m telescopes 
deliver high SNR observations of low surface brightness ob- 
jects and the NIR spectrographs available on ground based 
telescopes are better able to probe obscured dusty regions. 
We have undertaken a careful study to identify targets with 
bright reference stars close to the galaxy centre which are 
under represented in the current M.-cr plane. 

1.2 NGC 1399 

The first of these targets to be observed is the giant elliptical 
NGC 1399, which is the most luminous galaxy in Fornax and 
has a core- like s urface brightness profile characteristic of a 
cD type galaxy llKilleen fc BicknellllTgsl ISchombertlll986ll . 
Many dynamical s tudies in the visible have been under- 
taken in the past jBicknell. Bruce. Carter fc KilleenI [l989|: 
D'Onofrio et alJll995l: iFranx. Illing worth fc Heckman''l989l: 
Graham et alJl998l : lLongo et alJl9 94: Saglia ct al. 2000) al- 
though none have been made in the NIR or with the spatial 
resolution available with an AO system. 

The velocity dispersion of the spheroidal component 
ff» is calc ulated in different ways by diff'erent authors 
jFerrarese fc M crrit 2000: Gebhardt 2000). As highlghted 
bv T02. iFerrarese fc Merrltjl feoodh and subsequently FF05 
calculate the RM S dispersion with in a circular aperture of 
radius re/8 while ICebhardd tOOd) and T02 use the (lumi- 
nosity weighted) RMS dispersion within a slit aperture of 
length 2re. Whether or not the different definitions affect 
the measured slope of the relation, it is important to mea- 
sure the corresponding values for NGC 1399 so like can be 
compared with like. 

The effective radius of NGC 1399 is reported to be 
40 " bv[d ^ Vauco uleurs e t alj lll99lll . Using the longslit data 
of iGraham et alJ (1199811 to give a slit aperture of length 
2re centred on the galaxy (±40") and the photometry of 
iLauer et alJ i2005l) . we calculate the luminosity weighted 
RMS dispersion along the slit aperture to be 317±3 km s~^ 
(the quoted error is an estimate of random error only) . NGC 
1399 is therefore at the top of the T02 M.-cr plane, with 
three other galaxies that anchor the relation (M87, IC 1459 
and NGC 4649). The predicted BH mass for NGC 1399 is 
8.7 X 10*Mo so the Sol would be 0"34 (33pc); much smaller 
than even the best seeing at the best ob servatory sites. 

FoUowing IFerrarese fc Merrid ll2000l) . the RMS disper- 
sion within a circular apert ure of radius re 78 was estimated 
from the central 2" x 5" of iGraham et al.l il998ll 's data to 
be 329 ± 4km s~^ (as before, the quoted error is a measure 
of the random error). The relation of FF05 then predicts a 
BH mass of 1.89 x 10^ and the Sol would be 0'.'78 (75pc): 
still not well resolvable with the best seeing conditions. 

Fortunately, NGC 1399 is ideally suited to AO assisted 
observations. A bright (m„ = 13.8) reference star exists 
only 17" 6 away from the galaxy centre. The galaxy is al- 
most spherical with little or no rotation, so the slit is free to 
be aligned to the galaxy centre and the AO reference star. 
This allows us to monitor the AO correction as a function 
of time. 

One problem that plagues all stellar-dynamical esti- 
mates of black hole masses is the degenerac y between 
mass and anisotropy jBinnev fc MamonI lijill . Consider 
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the case of a spherical galaxy. Using the Jeans equa- 
tion, the mass enclosed withi n radius r can be written as 
iKormendv Richstonf^lToQ^h 



M{r) = -c+^ 
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(3) 



where v is the rotation velocity, cr^, ag, are the radial 
and azimuthal components of the velocity dispersion and j 
is the deprojected luminosity density. The first two terms in 
square brackets can be estimated almost directly from obser- 
vations. Both are positive for the vast majority of galaxies. 
The last two terms, however, depend on the the unknown 
anisotropy and can take either sign. Their effect on AI{r) 
is minimised for galaxies with steep j{r) profiles, steep ve- 
locity dispersion profiles ar and rapid rotation v ^ 0, all 
of which tend to be satisfied in power-law galaxies. Core 
galaxies like NGC f399, however, tend to be non-rotating 
with shallow density and dispersion profiles. For such galax- 
ies it is particularly important to constrain the anisotropy 
by modelling at the detailed shape of the galaxy's line-of- 
sight velocity profile s (VPs), f or wh ich high signal-to-noise 
spectra are essential iGerhardI lll993f) . 

Using the NAOS AO system coupled with the CON- 
ICA NIR imager / spectrograph at the European Southern 
Observatory's Very Large Telescope (ESO VLT) we have 
resolved the Sol of NGC 1399 and measured its stellar kine- 
matics using the CO absorption bands at 2.3/jm and the Cal 
absorption feature at 2.26/im. The SNR of the spectra range 
from ~ 90 to ~ 20 over the region used to extract kinemat- 
ics, with a SNR of ~ 70 at the CO bandhead (2.3/im). We 
use these kinematics to construct a spherical orbit super- 
position model for the galaxy to estimate the mass of the 
central MDO. Throughout this paper we assu me a distance 
of 19.9Mpc to NGC 1399 jTonrv et al.ll200ll) : the reader is 
reminded that BH mass scales linearly with assumed dis- 
tance. 

The structure of this paper is as follows. The data re- 
duction techniques are discussed in Section |21 the kinematic 
analysis is discussed in Section |3 the imaging and kinemat- 
ics are presented in Section 3] and the discussion of their 
implications is contained in Section |K| The dynamical mod- 
elling is described and discussed in Section|S| Finally, Section 
Q concludes. 



2 DATA AND REDUCTION 
2.1 Observations 

AO assisted K band images (Ks filter) and K band long- 
slit spectra (SK filter) of the nuclear region of NGC 1399 
were obtained wit h NAOS-CONICA jRousset et al.lll998l: 
iLenzen et ai]|l998l) on the nights of 30/11/03 and 01/12/03. 
The spectral range extended from 1.79/.im-2.45^m, with a 
scale of 0.972nm per pixel, although atmospheric transmis- 
sion limits high SNR data to 1.95/im-2.45/.im. The spatial 
scale of the spectroscopy was 0'.'0543 per pixel; the scale of 
the imaging was 0'.'027 per pixel. The slit width was 0'.'172 
(17pc) corresponding to an instrumental resolution (A/AA) 



10'° 
108 
10= 
10' 
108 



* stars 

■ ionized gas 



# y 



. . . . I 



_L 



' ' ' ' I 

o ellipticals 

♦ lenliculars 

? spirals 



V 



50 100 500 50 100 500 

a, (km s"') (7. (km s~') 

Figure 1. The sample of T02 and their best fit correlation of 
the form of Q with a = 8.13 and /3 = 4.02. Left: the symbols 
indicate the technique used to derive the black hole mass. Right: 
the symbols represent the different morphological types. 



of 880 at 2.3fim (as measured from the width of the arc lines) 
and an instrumental broadening of atnst = 145km s~^. The 
conditions while observing were excellent with the seeing 
varying between 0'.'4 and O'.'G. 

A total of 27000 seconds (7.5 hours) of useful on-source 
spectroscopic integration was achieved with individual ex- 
posures lasting 300 seconds each. Airmass ranged from 1.0 
to 1.4, although 85% of exposures were made with airmass 
< 1.2. The standard ABBA technique of nodding back and 
forth along the slit removed the need for separate sky expo- 
sures and Fowler readout mode was used to minimise read- 
out noise. AO assisted images of NGC 1399 were also taken 
for a total of 80 seconds with the Ks filter together with an 
equivalent number of sky exposures. 

In order to remove the complex atmospheric transmis- 
sion curve (referred to as telluric absorption) from the object 
spectra, it was necessary to observe several telluric standard 
stars (see Sec. 12.31 which have almost intrinsically feature- 
less spectra in the NIR (e.g. very hot O or B stars). As 
the strength of telluric absorption depends on the airmass, 
on the first night telluric standards were observed at air- 
masses of 1.02 (HD25631) and 1.40 (HD25631), and on the 
second night telluric standards were observed at airmasses 
of 1.12 (HD480) and 1.25 (HD41814). The difference in air- 
mass between telluric and galaxy exposures on each night 
was never more than 0.2 airmasses. Two kinematic template 
stars were also observed at low airmass; HD11931 (K4III) 
and HD25840 (MOIII). 

All stars (telluric standards and kinematic templates) 
were observed with an identical spectrograph configuration 
to the galaxy observations. However, for the stellar observa- 
tions alone, it was necessary to reduce the higher order gain 
of AO correction to ensure that the FWHM of the point 
spread function (PSF) was larger than the width of the slit 
to match the spectral resolution of the galaxy and stellar 
observations. With full AO correction, the FWHM of the 
PSF was ^O'.'l which would have been significantly smaller 
than the 0'.'172 wide slit. 

The position angle (PA) of the slit was 5.06° so as to 
include the galaxy centre and AO reference star in the slit. 
This allowed us to monitor and asses the AO correction. We 
assume that the SMBH lies at the most luminous region of 
the galaxy (also assumed to be at the centre of the galaxy) 
so it is important to position the slit to sample this region or 
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to be able to quantify any offset from it. Prior to acquiring 
the galaxy, an image of the slit on the detector was made 
by removing the grism and illuminating the slit with the 
flat field lamp. The slit image was then used as a bias for 
the subsequent acquisition images to precisely align the slit. 
Special care was taken with the first acquisition of the first 
night to ensure the PA of the slit intersected the AO refer- 
ence star and the brightest part of the galaxy. The PA was 
subsequently held fixed for all observations and only shifts 
perpendicular to the slit length (i.e. along the slit's minor 
axis) were made to maintain the slit position on the star and 
the galaxy centre. On the second night, the same PA was 
verified to hold the star and the galaxy centre in the slit and 
then shifts were made along the slit's minor axis as before. 



2.2 Reduction Techniques 

The data reduction was completed with the aid of the IRAF^ 
and ECLIPSE^ packa ges as well as custom IDL'^ scripts, 
incorporating use of the IDL Astronomy User's Library 
iLan dsman 1993)*. 

The spectroscopic data were initially reduced follow- 
ing the standard ABBA technique for NIR data reduction'' 
which has many advantages: the timescale on which the 
background subtraction is achieved is as short as possible, 
helping to correct for the variability of the NIR sky; any 
residual sky in a single A-B frame cancels with the residual 
in the B-A frame, assuming a uniform sky field; the pixel- 
to-pixel subtraction is very well suited to removing system- 
atic errors. However, such a technique does not optimise the 
random noise (sky, thermal, dark, readout) as a background 
exposure of the same duration as the object exposure is sub- 
tracted from each pixel. 

In an effort to increase the SNR, the background level 
of each pixel was interpolated as a function of time from 
all the data frames. Pixels with significant source (galaxy) 
counts were excluded when fitting a 3rd order polynomial as 
a function of time to each pixel position. In order to gauge 
the change in the SNR from interpolation, the random noise 
of the sky dominated region between galaxy and the AO ref- 
erence star was measured as a function of wavelength and 
compared to the noise of the frames without background 
interpolation: the interpolated background showed a signif- 
icant decrease in random noise (~ \/2 lower). This led to a 
significant increase in the final SNR of the data while main- 
taining the pixel-to-pixel subtraction to reduce systematics. 
However, the residuals from one object-sky pair do not nec- 
essarily cancel with the next pair, unlike the classical reduc- 
tion technique for ABBA sequence observations (assuming a 
uniform illumination). The interpolation is also susceptible 



^ IRAF is distributed by National Optical Astronomy Observa- 
tories (NOAO), http://iraf.noao.edu 

^ ECLIPSE is a reduction package developed 
by the European Southern Observatories (ESO), 
http: / /www. eso.org/projects/aot/eclipse/ 

Interactive Data Language, Research Systems, Inc. 
^ http://idlastro.gsfc.nasa.gov/ 

^ For further information on NIR data reduction, the reader is 
referred to the ISAAC data reduction guide (the NACO data 
reduction guide is still being constructed) 



to bias from bad pixels which can propagate into neighbour- 
ing frames. To counter this problem, the frames were cleaned 
of bad pixels prior to (object and sky frames), during (sky 
pixels only) and after the interpolation (background sub- 
tracted pixels). Indeed, a consequence of interpolating the 
background was the identification of faulty pixels on the de- 
tector and the construction of an very accurate bad pixel 
map. 

The timing and conditions of the spectroscopic data 
make it well suited to sky interpolation: 7.5 hours of obser- 
vations were performed on neighbouring nights with limited 
interruption between the exposures on each night and the 
atmospheric conditions during the two nights were excellent 
and stable. 

After correction for the odd-even effect of the detec- 
tor, sky subtraction, flat fielding, bad pixel correction, field 
transformation and wavelength calibration, the exposures 
were aligned along the spatial axis (to an accuracy of a few 
tenths of a pixel) by using the centroid of a Gaussian fit- 
ted to the light profile of the AO reference star. Reduction 
of telluric standard star spectra followed the same sequence 
except with only standard ABBA sky subtraction. 

The image data was reduced in the standard manner 
for NIR observations. After correcting for odd-even effect, 
the sky exposures were subtracted from the object frames, 
followed by flat flelding. The centroid of a Gaussian, fitted to 
the AO reference star, was then used to align the individual 
exposures to a tenth of a pixel accuracy. 



2.3 Telluric Correction 

Molecular gas in the Earth's atmosphere (CO2, H2O etc.) 
produces an absorption spectrum known as telluric absorp- 
tion and this must be removed from astronomical obser- 
vations. The galaxy and kinematic template observations 
were corrected for telluric absorption in the standard man- 
ner iOriglia. Moorwood and Olival Il99'3) . the details and 
caveats of which are described below. 

The stellar continuum must first be removed from the 
telluric standard star spectrum. We therefore divide through 
by a blackbody spectrum of appropriate temperature to the 
star. This is a blind division; one cannot fit the temperature 
as the atmospheric transmission curve is still imprinted on 
the spectrum (in addition to any instrumental response). 

Telluric absorption varies with airmass so the strength 
of the sharp prominent absorption features in the normalised 
telluric spectrum must be matched to those in the object 
spectrum (of similar airmass); division of this optimised 
standard then removes the prominent (high frequency) fea- 
tures and, we assume, the low frequency component of the 
atmospheric transmission curve. The telluric spectrum, now 
free from stellar continuum, is not yet normalised so we di- 
vide by a linear flt to regions where we believe the atmo- 
spheric transmission to be ~100%: (2.035-2.04)/im, (2.09- 
2.15)/xmand (2.21-2.22)/im. The increased strength and oc- 
curance of telluric features beyond 2.3/im causes atmo- 
spheric transmission to be consistently below unity, and 
falling with increasing wavelength so one must extrapolate 
the continuum from shorter wavelengths. The absorption 
depth in the now normalised telluric spectrum, A{X) is then 
amended using the expression 
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A'{\) ^1.0- F[l.O- A{X)] (4) 



where F is a free parameter, to account for small variations 
in airmass between the object and telluric standard star ob- 
servations. Note that this linear correction for airmass is 
only an approximation to a more complicated response and 
is only effective for very small differences in airmass. The ini- 
tial normalisation was then removed by multiplying back by 
the original linear fit. Every telluric spectrum also required a 
small shift in wavelength (around a few tenths of a pixel) to 
compensate for differences in wavelength zero-point between 
the object and telluric spectra. The optimisation of airmass 
and wavelength was initially automated by minimising resid- 
uals around the strong telluric features at (2.0 - 2.1)/xm in 
the telluric divided object spectrum. Fine adjustments were 
then made to further optimise the removal of prominent tel- 
luric absorption features at the region of interest (the CO 
band-heads after 2.3/xm). Note that in practice, while one 
can correct individual stellar exposures (due to the high SNR 
in each frame) , galaxy exposures must be coadded in groups 
of similar airmass to increase the SNR before an accurate 
telluric correction can be determined. 

Telluric correction is not without complication though. 
The telluric correction is optimised for the removal of sharp, 
prominent, high frequency absorption features. Accordingly, 
there is likely to be minimal residuals from such in the cor- 
rected spectra. However, variation in the continuum normal- 
isation of a galaxy spectrum is known to introduce system- 
atic e ffects into the derived kinematics jvan der Marel et alJ 
Il994h . Due to the blind division of a blackbody spectrum and 
the lack of continuum reference in the telluric spectra after 
2.3fim, error in telluric correction will most likely manifest 
itself as error in the continuum normalisation of the telluric 
spectrum after the application of Q • This would propagate 
into an error in the continuum level of the object spectrum 
(galaxy or kinematic template). Such an error may not be 
uniform over the length of the spectrum: there is likely to be 
a higher chance of error where the normalisation was com- 
pletely extrapolated, after 2.3/im. To help compensate for 
such systematic effects, we include a polynomial continuum 
correction when extracting kinematics (see Sec. |3Jl, but the 
constraint for the correction is the minimisation of Xs 0) 
which does not guarantee to choose a solution free of sys- 
tematics. Furthermore, this continuum correction is additive 
whereas any real error would be divided into the object spec- 
trum. 

Different functions were used to normalise the galaxy 
and kinematic template spectra after telluric correction. The 
galaxy continuum appeared linear with wavelength and so 
at every position along the slit, a linear fit to the continuum 
shortward of the CO bands was sufficient. The continuum 
of each kinematic template shortward of the CO bandhead 
was clearly non-linear but was well fitted (and removed) by a 
blackbody spectrum. As the continuum after 2.3/im must be 
extrapolated, it was necessary to fit slowly varying functions 
with relatively little freedom. Note that with all continuum 
fits, care was taken not to include obvious absorption or 
emission features, or areas of significant telluric absorption. 



2.4 Kinematic Templates 

The choice of kinematic template(s) is important to accu- 
rately extract kinematics from the galaxy spectra. The sys- 
tematic errors introduced into kinematics by use of a poor 
template are well studied (van der Marel et all If 993) . al- 
though it is difficult to numerically quantify such effects. 
At best, we can say that different templates appear to in- 
troduce systematic offsets into the VP parameters. Ideally, 
a large library of spectral types should be available so that 
an optimal mix can be found, but in this case, due to time 
constraints, only two templates were observed with the same 
instrumental setup as the NGC f 399 observations. However, 
it is possible to check if the templates are well matched to 
the luminosity weighted population of the galaxy. 

The CO (2-0) band head at 2.2935/im is an indica- 
tor of stellar type: a linear relation has been found be- 
tween the equivalent width (EW ) of the first CO feature , 
Wgnf^-m and the stellar type l|Kleinmann fc Halll If 9861 : 
lOriglia. Moorwood and Oliva '' f99jl. Although this relation 
is based on observations of individual stars, it can be ex- 
tended to galaxy populations by applying a correction based 
on the velocity dispersion of the system jOliva et alJlf995l : 
Thatte, Tccza and Gcnzol 200(|). The relationship differs be- 
tween giants and supergiants but this should not be a prob- 
lem for the old, giant dominated, population of an elliptical 
galaxy such as NGC f399. To estimate the galaxy EW, the 
(rest frame) wavelength range over which the CO (2-0) EW 
is defined must be shifted to the velocity frame of the galaxy 
and the EW measurement must be corrected for the galaxy's 
velocity dispersion. The kinematic properties of the galaxy 
are not known a priori and depend on the template used to 
extract them. However, it is possible to account for this and 
calculate reasonable limits on the EW of the galaxy. 

Such analysis was performed for NGC f 399. The vari- 
ation of EW with velocity dispersion was simulated using 
the kinematic templates and the necessary quadratic cor- 
rection (in a) found; this was then applied to the galaxy 
EW measurements. Furthermore, all measure ments were 
scaled by a constant factor (the correction of lOliva et alJ 
if 995ll l to correct for the instrumental resolution of CON- 
ICA, allowing direct comparison with the relation of 
lOriglia. Moorwood and Clival (If 9931) . The results are pre- 
sented in Fig. H] and discussed in Sec.|ll 

2.5 AO correction and PSFs 

The quality of the AO correction can be estimated from the 
reference star, which was observed in the slit simultaneously 
with the galaxy. The f D profile of the star (calculated from 
summing the fiux over the same wavelength range that the 
kinematics are extracted from) is well fit by a double Gaus- 
sian as shown in Fig. (21 

However, the correction and PSF vary further away 
from the reference star. Although the exact 'off-source' PSF 
at the centre of the galaxy is unknown, it can be estimated. 
Using the NAOS preparation software vf .74^, one is able to 
simulate how the PSF varies with seeing, airmass and dis- 
tance from the reference star. However, the simulated PSFs 

^ http:/ /www.eso.org/observing/etc/naosps/doc/NAOS-PS- 
tool.html 
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1.214 


1.280 


72 


0.010 


0.129 


0'x2 


3.667 


3.501 




2.784 


2.630 


FWHM along x (") 


0.135 


0.150 


FWHM along y (") 


0.181 


0.189 


Strehl Ratio (%) 


38 


30 



and off-source PSFs: 
T<^x2,o-y2) of ^ 2D dou- 



Table 1. The characteristics of the 
the best-fit parameters (71 , cr^,!, o-yi, 72, 

ble Gaussian jSJ fit to the PSFs after incorportaing broadening 
effects (illustrated in Fig|2]and subsequently used in dynamical 
modelling); the FWHM of both broadened simulated PSFs along 
the X and y axes; the Strehl ratios of the simulated PSFs prior to 
incorporating additional broadening. 



do not account for slit effects, non-perfect data reduction 
(such as error in the field transformation or frame ahgn- 
ment) or drift of the the tip-tilt correction during long ex- 
posures. AU these effects will further broaden the PSF. To 
account for such, we convolve simulated 2D on-source PSFs 
(with various airmass and seeing conditions) with a, 6 x 3.2 
top hat (3.2 pixels is the slit width and 9 accounts for im- 
perfect data reduction and tip-tilt drift). We optimise 6 and 
the atmospheric conditions by minimising the difference be- 
tween the ID profile along y = and the observed ID profile 
of the AO reference star (Fig.l^J. Good agreement between 
the simulated and observed PSFs is found for 0'.'6 seeing, an 
airmass of 1.2 and 9 = 2.2 pixels. To estimate the off-source 
PSF at the galaxy centre, we simulate the 2D off-source PSF 
with identical atmospheric conditions and convolve it with 
the same 2D kernel (we assume that the broadening effects 
are uniform over the field of fiew) . The resulting on- and off- 
source PSFs are both well fitted by a double 2D Gaussian 



PSF{x,y) = 



+ 



exp 



"2 (<T,i) + (^/i) 



■ exp 



+ 



y 

CT„2 



(5) 



The best-fit parameters are given in Table Q and are used 
to describe the 2D PSF in the dynamical modelling. The 
Strehl ratios also quoted in Tableware calculated from the 
unbroadened (pre-convolution) simulated PSFs: it makes 
no sense to compare the ideal Airy disk pattern with the 
slit-convolved PSF of the spectrograph. The FWHM of the 
off-source PSF along the x-axis is 0'.'15 (2.75 pixels), corre- 
sponding to 14pc, which we adopt as our formal resolution. 

The effect of contamination from the seeing-limited halo 
on the spectra and the knock-on effect on the VPs is not 
known precisely. According to Q, if the galaxy spectrum 
is a weighted sum of many different spectra with many dif- 
ferent VPs, the problem is linear and Xs be minimised 
for a similarly weighted mean VP. The galaxy spectra and 
kinematics are therefore only expected to be 'diluted' by this 
effect. However, the dynamical modelling uses the estimate 
of the off-source PSF to fully account for the shape of the 
PSF; the peculiar shape of the AO corrected PSF does not 
bias the models or the derived BH mass. 

The quality of the PSF was assessed as a function of 
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Residual 




0.0 

Radius (arcsec) 



Figure 2. The profile of the AO reference star together with a 
double Gaussian fit and the residuals from this fit. The best-fit a 
for each of the Gaussians are 0'.'048 and 0'.'122 and the total fluxes 
(7s) are in the ratio 2:1 respectively. This ID slit profile of the 
AO reference star is used to estimate the 2D on- and off-source 
PSFs. 



time from the AO reference star in the slit. The quality of the 
overall AO correction for the final coadded galaxy data can 
be marginally improved with frame selection. However, the 
cost in SNR was too high for any significant improvement 
in this case, so frame selection was not implemented. 



3 KINEMATICS 

We assume that each galaxy spectrum G(A) is a Doppler- 
broadened version of some underlying "average" stellar tem- 
plate T{X). If we could completely remove all continuum 
features from G and T, then we would have that 



G{X) ^T{X)(g) L{v), 



(6) 



where L{v) is the unknown line-of-sight stellar velocity pro- 
file and (8) denotes convolution. In reality, it is impossible 
to completely remove all continuum contamination from ob- 
served spectra by eye, nor do we know the correct template T 
to use. Furthermore, both G and T are measured with finite 
signal-to-noise, making it impossible to use a simple decon- 
volution method to obtain L{v) directly from 0. 

Many methods have been pr oposed to address these 
proble ms. We follow 'Rix fc White! lll992h . lSaha fc William^ 
l|l994r) and Ivan der Marel (11 1994) and extract L{v) using a 
direct pixel-fitting method. Taking the continuum-divided 
star and galaxy spectra we find L{v) and its associated un- 
certainties by minimising 



E 



AG. 



(7) 



where AGi is the measurement error in the ith galaxy spec- 
trum pixel, A^p is the number of pixels being fitted and the 
third term corrects for any continuum that escapes our ini- 
tial continuum division. The parameter fc ~ 1 accounts for 
differences in normalisation between G and T: we want our 
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VPs to be normalised with J L{v)dv = 1. We lurther as- 
sume that T(A) is well approximated by a weighted aver- 
age of known stellar templates. As we only have two stellar 
templates available for kinematic extraction, we define our 
optimal template to be 



T(A) = /Ti(A) + (l-/)r2(A) 



(8) 



where / defines the relative fraction of each of the two avail- 
able templates (Ti, T2). 



Richstonc & SarEent' 
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'Franx. lllineworth & HeckmanI 


1989 




Bender 


1990; 


van dcr Marcl & Fraiix 


„199al, pixel 



fitting does not require any assumptions about window 
functions and allows us to propagate the pixel-to-pixel 
error estimates in our measured spectra directly into 
uncertainties in our LOSVDs. Note that the Nal doublet 
at 2.21^m does not appear to be well fit by our templates 
(the line strength in the galaxy spectra is much higher 
than that in our templates). Thus the fitting range for 
extraction of kinematics is from 2.249/im to 2.438/xm (rest 
frame) to include the Cal feature and the CO bands. We 
use two different parametrisations for L{v), each of which 
is discussed separately below. 



of the expansion ^ from the VP L{v). Since the Hi form a 
complete set, we can make Xo arbitrarily small for smooth 
L{v) simply by increasing the number of terms A'^ included 
in the series, although the expansion © is of course useful 
in practice only if it can provide a good fit for small N. 
We note that when minimising Xo fixed (7, V, a): 

(i) the Gauss-Hermite coefficients hi are modified mo- 
ments Ijllfl of L(v); 

(ii) the Hessian d^Xo/9hidhj is diagonal, meaning that 
the hi are independent; 

(iii) the hi are also independent of the choice of N; 

(iv) if we choose (7, V, a) to be the parameters of the 
Gaussian that minimises 1121 then {ho,hi,h2) = (1,0,0) 
and (/i3,/i4) measure the lowest order anti-symmetric and 
symmetric deviations from this Gaussian. 

Because of this last point, most Gauss-Hermite parametri- 
sations of VPs fix {ho,hi,h2) — (1,0,0) and use 
(7, V, a, /13, hi, . . . , /iat) as the free parameters in the fit. For 
our purposes, this approach has the drawback of introduc- 
ing nontrivial correlations between the Gaussian parameters 
(7, V, a) and the {hi}, making these parameters cumbersome 
to use when comparing to dynamical models. 



3.1 Gauss-Hermite parametrisation of the VP 



3.1.2 Real data 



VPs are expected to be reasonably close to Gaussian. 
A convenient way of parametrising a VP L(v) is by us- 
ing a truncated Gauss-Herm ite expansion iGerhardI Il993l : 
Ivan der Marel fc Fran^ll993l) . 



Lgh{v) 



2-Ka 



exp 



in which one starts from a Gaussian with scale factor 7, 
mean V and dispersion a and uses a weighted sum of Her- 
mite polynomials Hi to quantify deviations of L{v) from 
this underlying Gaussian. This parametrisation is degener- 
ate: there is a different set of coefficients {hi} for each choice 
of (7, V, a) and a number of different procedures have been 
used to fit spectra using @. Before describing our own pro- 
cedure, we review the motivation for these earlier methods. 



3.1.1 Perfect data 

Consider first an idealised situation in which we knew L{v) 
perfectly. Then there is a unique set of expansion coefficients 
{hi} for any (sensible) choice of Gaussian (7,!/, cr). To see 
this, recall that the Hi satis fy the orthogonality relation 
Jvan der Marel fc Fran^ll993ll . 



1 / 1 

— / exp{-x')Hi{x)H.j{x)dx = —^Sij, 

from which it is easy to see that choosing 
1 



hi = 



\/27, 



poo 




/ exp 


2 V o- J 


/ —00 





'v-V 



minimises the mean-square deviation 

/> oc 

' [L{v) ^ LG-a{v)f dv 



xl 



(10) 



du(ll) 



(12) 



In reality we do not have direct access to L{v). Instead 
we constrain it by investigating how well a parametrised 
form, such as 0, affects the fit to the discretely sam- 
pled gala xy spectrum Gj usi ng the Xs given by JTJ. As 
^an dcr Marel fc Franxl ^1993^ point out, eq. reduces 
to 1121 1 in the limit of high resolution, finely sampled spec- 
tra and sharp template features. Here we consider the case 
where this limit does not apply. Let us assume that the con- 
tinuum has been perfectly removed from d and that T is 
the correct stellar template. Then ||7J becomes 



Xs 



E 



{T ® {L ~ LG-a))i + Ui 



AG, 



(13) 



where L{v) is the galaxy's real underlying VP and Ui = 
Gi — {T ® L)i is the noise in the i**" pixel. Now if we fix 
(7, V,a) and minimise 1131 with respect to the {hi}, then: 

(i) in the absence of noise, the coefficients hi axe still the 
modified moments 1111 of L{v); 

(ii) the hi are not independent since the Hessian 
d^Xa /dhidhj is no longer diagonal; 

(iii) since the hi are not independent, in the presence of 
noise there is a different set of hi for each choice of N; 

(iv) if we choose (7, V, a) to be the parameters of the best- 
fit Gaussian to (O, then the minimum Xs will not occur at 
precisely (/iq, /ii, /12) — (1,0,0). 

Our procedure for fitting Gauss-Hermite coefficients is 
motivated by these points and by our desire to have a set of 
parameters that depend linearly on L(v). The procedure is 
as follows: 

(i) Choose (7, V, a) to be the parameters of the best-fit 
Gaussian to the VP: find (7,^,0"), template fraction / and 
continuum parameters c; that minimise ||7| with k — ho = 1 
and hi = h2 — ■ ■ ■ — 0. 
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(ii) Holding (7, V, a) fixed at their best-fit values, we find 
the hi, ci and / that minimise l|7J. Having found these hi, 
the normalisation is fc = 7(^0 + h-^l^pl + /i4-y^3/8 + • ■ •). 

(iii) Finally, we divide 7 and the hi by fc. 

We use a standard Levenberg-Marquardt routine to 
carry out the minimisations in the first two steps. Our best- 
fit parameters are (fto, ■ • • , ^iv), along with their covariances 
and the choice of (7, V, ct). Note that there are no errors as- 
sociated with (7, V, cr) in our version of the Gauss-Hermite 
parametrisation: they merely reflect the Gaussian around 
which we have chosen to expand L{v). Choosing the best-fit 
Gaussian here lets us make a straightforward comparison of 
our kinematics with earlier work. In practice we find that 
our method yields (/lo, fei, /12) that differ from (1,0,0) by 
around (0.1,0.04,0.04), but with very strongly coupled er- 
rors among all the even hi. We describe how we deal with 
these covariances in section below. 

We have tried using standard simultaneous A'^ -|- 1 pa- 
rame ter fits to (7, y, (J, /13, . . . , /ijv) ivan der Marel fc Franxl 
Il993l) . but our VPs are so strongly non-Gaussian at the cen- 
tre of the galaxy that the usual linea r approximations among 
the errors in these parameters (Ivan der Marel fc Franxl 
Il993l) break down and we find multiple minima in Xs 0- 
For spectra at (-0'.'08,0'.'02) this process can yield /14 > 0.6 
and alarmingly low values of a (~ 250km s~^) which de- 
scribe a triple peaked VP. In fact, for such values of /14 it 
can be shown that even the idealised Xo II12I I has multiple 
minima in a and /14. 

3.2 VP histograms 

One might expect that the putative BH in NGC 1399 
would cause high-velocity wings in the central VPs, which 
might not be captured well by the low-order Gauss- 
Hermite parametrisation above. Therefore we also fit "non- 
parametric" VPs, in which we choose A'^ regularly spaced 
velocities v\ < v-i < ■ ■ ■ < vn and parametrise L(v) as the 
histogram 

Ly^[v)=^US^{v), (14) 

i = l 

where the step function Si(v) = \ \i Vi < v < Ui+i and is 
zero otherwise. 

Given parameters L\,. .. , L„„ it is straightforward to 
calculate the convolution Q of this Lh(w) with a stellar tem- 
plate. For any given galaxy spectrum G there will be many 
sets of parameters that produce good fits to the spectrum, 
but most of them will be unrealistically jagged. Therefore, 
instead of minimising the Xs given by eq. directly, we 
minimise the penalised Xp = Xs + P\^i\^ where the penalty 
function 

P[L,] = a^(L.+i-2L, +L._i)' (15) 

i 

uses the mean-square second derivative of L{v) as a measure 
of the jaggedness of the solution. 

Our procedure for fitting Li is simple: 

(i) Find the best-fit smooth Li, continuum parameters c; 
and template fraction / by minimising the penalised Xp = 
Xs + P\U\ with fe = 1. 



(ii) Set the normalisation factor to fc ^ = ^^(n^+i — iii)Li 
and rescale the Li. 

This makes no attempt to impose the obvious non- 
negativity constraint on the Li. While the resulting VP his- 
tograms are fine for "by-eye" comparisons of one VP against 
another, they are not suitable for direct comparison against 
models. So, in >l6.3l we describe a variation on this fitting pro- 
cedure that takes account of the correlations among the Li 
and removes the bias introduced by the penalty function. 

Based on an average separation of 300A between the 
CO bands in the fitting range (1^00(2 - 0), ^^CO(3 - 1), 
^^CO(4 - 2), 1^00(5 - 3), CO (6 - 4)), the maximum 
relative velocity we can reasonably hope to measure is 
~ 1900km s-\ The systematic velocity of NGC 1399 is 
~ 1500km s~^ so we divide each LOSVD into 71 = 50 eq- 
uispaced velocity points between v\ = —1000km s~^and 
Vn — 4000km s"'^. We choose a = 4 x 10^, which is the min- 
imum value required to give smooth non-parametric VPs 
consistent with our outer Gauss-Hermite VPs. 



4 RESULTS: IMAGING AND KINEMATICS 

An 80 second, AO corrected, Ks exposure of NGC1399 is 
shown in Fig. |3 It has the same PA as the long-slit ob- 
servations and a slit image (0'.'172 wide and centred on the 
AO reference star) has been overlaid. If the reference star 
is perfectly centred in the slit, there is a maximum posi- 
tion error of 0.3 pixels (1.6pc) on the brightest region of 
the galaxy from error in the PA. However, the acquisition 
images (incorporating flat fielding etc. and using the slit im- 
ages) indicate that the centring of the star along the minor 
axis of the slit was accurate to only a pixel (5.2pc); the slit 
was approximately 3 pixels wide (0'.'172 or 17pc). Hence, the 
dominant error in aligning the slit on the galaxy centre is 
from shifts along the minor axis, not from error in the PA. 
The former is expected to be random about the centre of 
the slit; the later would give a systematic difference for all 
observations. The final position accuracy is therefore good: 
the slit was aligned on the brightest point of the galaxy to 
an accuracy of around ±0'.'06 (±6pc), with no significant 
systematic offset. 

Two globular clusters are also seen in the image data 
close to the centre of the galaxy. The nearest is 1'.'15 to the 
east of the galaxy centre. A Gaussian fit to this globular 
cluster yields an estimate on the FWHM of the image PSF 
to be 0'.'078 (7.5pc). This is equal to the on-source FWHM 
predicted by the NAOS preparation software in Sec. 12.51 
suggesting that the single Gaussian was probably fit to the 
diffraction limited component of the globular cluster pro- 
file and ignored the seeing limited halo, although the atmo- 
spheric conditions were excellent around the time of obser- 
vation, with the seeing occasionally dropping below 0'.'4. 

The parametrised kinematics of the central few arcsec- 
onds of NGC 1399 are sho wn in Fig. |H Where possible, pre- 
vious data is over plotted dCraham et aljll998t iLongo et alJ 
[T994I . The fiux is calculated by summing over all available 
wavelengths at each point along the slit. The parameters 
(7,w,(j) are chosen to be the best-fit Gaussian parameters 
and the Gauss-Hermite coefficients hn are derived using 
these best-fit Gaussian parameters. The velocity v is given 
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Figure 3. Top: A Ks band image of NGC 1399 showing the slit alignment (white) on the galaxy and the AO reference star (saturated 
on this colormap); the PA of the slit is 5.06° so that the reference star is approximately due north of the galaxy nucleus; the nucleus of 
the galaxy and the star are separated by 17''5. Bottom: the same as the above but magnified and centred on the nucleus of NGC 1399 
with isophote ellipses over plotted in blue. Note the elongation of the nucleus to the south-east. Each pixel is 27mas wide (the pixel size 
of the spectroscopic data was twice this at 54maa). The angular scale is given in arcseconds. 



relative to the mean heliocentric velocity of 1467 ± 4 km s^^ . 
Note the error quoted for this velocity is an estimate of the 
random error only. 

Fig. 1^ compares the non-parametric and parametric 
VPs for the central arcsecond of the galaxy. Error bars are 
given for the non-parametric VPs. The galaxy spectra and 
best-fit broadened template (constructed from the Gauss- 
Hermite VP) for the central arcsecond are also shown. 



4.1 Highlights 

Many interesting features are present in the above figures. 
The parametrised kinematics show a decoupled velocity 
structure across the galaxy centre, a double peaked veloc- 
ity dispersion across the centre separated by 0'.'2 or 19pc 
(with the dip in a half a pixel off the galaxy centre) and 
significant variations in the hi with radius. The imaging 
shows offset asymmetric isophotes at the galaxy centre and 
the non-parametric VPs (while generally in good agreement 
with the parameteric VPs) show strong high velocity wings 



at r = — O'.'OS and lopsided velocity structure at r = 0'.'14, 
which are not detected with the parametrised fit. 



4.2 Further Details 

The kinematic features described above persisted when the 
spectra were reduced using the standard sky subtraction 
technique for ABBA sequence observations (rather than us- 
ing an interpolated sky), albeit with increased noise. Angu- 
lar distances in Figs.|31and|S]have been shifted such that the 
brightest region of the galaxy lies at r = 0". This position 
was det ermined by scaling to the (outer) centroid of a Nuker 
profile llLauer et alJll995li fitted to the fiux of Fig. H The 
photometry is not symmetric at the centre which could bias 
the calculation, although when the inner 0'.'5 of the light pro- 
file were omitted the change in the centroid was less than 0.1 
spectroscopic pixels (5mas or 0.5pc). The centroids of the el- 
lipses fitted to Fig.|3|are not constant and fluctuations of up 
to 0.5 image pixels (13.5 mas or 1.3pc) from the mean are 
present. The scale on Fig. |3|is aligned to the mean x-axis 
centroid of the ellipses. Registration between the angular 
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Figure 4. The nuclear kinematics of NGC 1399: The top plot 
indicates the fiux received at each position on the slit, the next 8 
plots show the variations in the parameters of |2j and the lower 
two show the template and galaxy EWs and the optimal kine- 
matic template composition as a function of radius from the nu- 
cleus ^_Wiiere_a£plicable, previous data is s hown as orange trian - 
gles iLongo et alj|l994l) and blue squares iGraham et al.lll99Sll . 
In the EW plot, the kinematic template EWs are shown in orange 
(MOIII) and blue (K4III); the width represents a 20km s~^ error 
in V and crinst- The EW of NGC 1399 is in shown in black (red 
shows the combined effect of 20km s~^ errors in v and a). The 
PA of the slit was 5.06° so that the positive radius is northward. 



See seperate figure 



Figure 5. The spectra and VPs of the central arcsecond of NGC 
1399. The centroid of the light profile along the slit (taken to 
be the nucleus) is roughly in between two pixels allowing for 
the VPs to be constructed at approximately equidistant inter- 
vals to the left and right of the nucleus, without interpolation. 
For every point along the slit, two VPs are plotted: black is the 
Gauss-Hermite VP jSJ and blue is the non-parametric VP. The 
velocity axis of the VPs has been corrected to the mean sys- 
tematic velocity of the system. Accompanying every spectrum is 
the best-fit broadened template (green) reconstructed from the 
Gauss-Hermite VP. 



scale on the image and that of the kinematics may not be 
perfect as the zero points were determined independently. 
However, is likely to be no worse than 0'.'05. 

EW measurements for the two kinematic templates and 
the galaxy are shown in Fig. |1] The template EWs (orange 
for the MOIII and blue for the K4III) allow for ±20km s"^ 
errors in systemic velocity and instrumental dispersion. The 
corrected EW of the galaxy is shown in black while the effect 
of ±20km s^^ errors in velocity and dispersion are shown 
in red. The asymmetry in the error is associated with the 
error in velocity: large positive or negative changes in veloc- 
ity move the CO (2-0) feature out of the predefi ned wave- 
length limits lOrid ia. Moorwood and Olival il993l) . thus re- 
ducing the EW. The kinematics {v, a) used to correct the 
galaxy EW are those shown in Fig. 01 Systematic differences 
in (i;, a) from using different templates in kinematic extrac- 
tion are less than 20km s^^ in each case. These EW mea- 
surements indicate that our kinematic templates are well 
matched to the luminosity weighted stellar population of 
the galaxy. With reference to lOriglia. Moorwood and Olival 
( 1993), one can see that the corrected EW measurements of 
the two templates (11.3 and 12.5 for the K4III and MOIII, 
respectively) are typical of their spectral classifications. The 
galaxy EW (average value of 12.1) is closer to that of the 
MOIII template; this is also reflected by the favouring of the 
MOIII template in the kinematic extraction (Fig. 2J ■ How- 
ever, one should be apprehensive about the variations of the 
galaxy's EW with radius: the dispersion correction assumes 
a Gaussian VP which we have seen is not always the case. 
Overall, the analysis of EWs indicates that the templates 
are well matched to the galaxy population, so we should 
not expect problems associated with template mismatch to 
be present in the kinematics. Additionally, previous (seeing 
limited) data on the central velocity dispersion of NGC 1399 
agrees with our dat a (Fig.^J, given the different resolutions. 
ISaglia et alJ i2000h report slightly negative /14 towards the 
galaxy centre which we confirm although our rise in /14 in 
the very central O'.'S is beyond the (seeing limited) spatial 
resolution of their data. 

To maintain a high SNR at larger radii (| r | > 0'.'3), 
the spectra were binned into pairs prior to kinematic extrac- 
tion. However, after binning, Xs rose, on average, from 
~160 to ~230 (an approximate change of \/2) for a fit to 
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174 data points. This indicates that systematic errors are 
starting to become comparable to the random errors. There 
are many possible causes of systematic errors. The optimal 
template, convolved with the best VP will not reproduce 
all the features in the galaxy spectrum; although the galaxy 
and template spectra have comparable CO (2-0) EWs, there 
may be absorption or even emission features in the galaxy 
spectra that are not obvious and are not accounted for in 
eq. Q, but still contribute to the differences between the 
galaxy and templates and to the systematic 'noise'. Telluric 
correction is only estimated to be accurate to 1%, increas- 
ing to 2%-3% at regions where the telluric absorption shows 
sharp prominent features. In addition, the stars used to ob- 
tain a telluric spectra are unlikely to have featureless spectra 
at high SNRs. The SNR of the data binned at C/.'S over the 
wavelength range used to extract kinematics falls from 110 
(2.25/im) to 30 (2.45/im) per pixel due to the thermal back- 
ground (the SNR at the unbinned galaxy centre falls from 
90 to 20). This SNR (and corresponding error estimates) 
quantify random error alone and do not account for system- 
atics introduced by telluric correction. Binning will reduce 
the random noise but not the systematic noise. Thus as Xs 
J?! is weighted by the random noise estimates, when the 
random error in each spectrum is reduced, it will increase 
so long as the systematic noise persists. 



5 DISCUSSION 

The interpretation of the results is discussed below in ap- 
propriate sections. 



5.1 Decoupled Kinematics 

There is a strong rotation gradient in v within a radius of 
0'.'5 (48pc) which is clearly decoupled from the kinematics at 
larger radii. The magnitude of the central rotation (taking 
the difference of the maximum velocity in each direction) is 
~70km s~^. The absence of high resolution data perpendic- 
ular to our long slit position prevents us from concluding 
if the system is truly counter rotating, or just decoupled. 
Previous pub lications considered the possibility of a d ecou- 
pled system dP'Onofrio et al.lll995l : ISadia et alJl2000h . but 
the poorer spatial resolution of the data prevented a reliable 
detection. Kinematically decoupled cores are often (but not 
exclusively) found in s pherical, high dispersi on systems with 
core- like photometry jEmsellem et al.(l2004^ . 



5.2 The Central Dispersion Profile 

The drop in a half a pixel off the galaxy centre is ac- 
companied by less convincing dips in 7 and the CO (2-0) 
EW. These features are close to the limit of the spatial 
resolution and thus may be unresolved. It is well known 
that random errors in 7 and a are statistically correlated 
jEfstathiou et al.lll980l) and these features are only clearly 
seen when continuum correction is incorporated into the 
minimisation of Xs 0i continuum correction is left out 
of the minimisation (i.e. calculated once from the contin- 
uum shortward of the CO (2-0) feature), 7 shows no ob- 
vious change (and is noisier in general) and the dispersion 
becomes flat at the centre. 



The decrease in 7 and a could be correlated effects from 
fitting constrained parametrised VPs to mo re complicated 
profil es; 7 is expected to fall at the centre jvan der MareJ 
119941) . However, the non-parametric VP at r = O'.'OS matches 
the parametric form well. Alternatively, the dip in a could 
reflect a genuine fall in the dispersion of the stars at the 
centre of the galaxy, but this is the exact opposite to what 
one would expect when approaching the BH. 

An alternative explanation is the presence of young stel- 
lar population at the centre, although 7 and a may not 
necessarily decrease in this case: a young stellar population 
may include younger MK types with reduced CO (2-0) EWs, 
but it would also have an increased fraction of supergiants, 
which have a highe r CO (2-0) EW than giants for the s ame 
MK spectral type dOridia. Moorwood and Clival Il993h . In 
addition, one would expect an correlated change in the op- 
timal template used to extract the kinematics, which is not 
the case here. One would expect a different stellar popula- 
tion to produce a colour gradient but unfortunately there is 
insufficient homogeneous data at the required resolution to 
confirm this. 

Central dips in a have been see n in other core galax- 
ies such as M87 a nd NGC 4649 (Ivan der Marel Il994l : 
iPinknev e^al]|2003^ . which both have nuclear activity. If 
NGC 1399 contained a cen tral nucleus of light from a non- 
thermal, non-stellar object (" van der 1^3X6111199411 . one would 
expect the CO bands to show a lower 7 relative to the con- 
tinuum level which may introduce correlate d dips in 7 and 
a. New photometry by iLauer et alJ l|200^ shows a slight 
excess of light at the galaxy centre after subtraction of the 
best-fit Nuker profile and therefore suggests the presence of 
a very faint nucleus. However, the authors emphasise the 
possible dangers and unknown systematics of extrapolating 
the Nuker fit into the central regions. 



5.3 Offset Photometry and Peculiar VPs 

The Ks band image of NGC 1399 shows a departure from 
spherical symmetry in the central 0'.'5 with an elongation 
of the surface brightness towards the east-south-east, ap- 
proximately 0'.'2 (19pc) in length. The image is a relatively 
short exposure and the SNR is low. However, the elonga- 
tion of the central isophotes covers more than 16 pixels; 
the probability of finding 16 neighbouring pixels above the 
mean, given the noise statistics of the image, is negligible 
so we can rule out the possibility of random error causing 
such an artifact. Obvious systematic errors which may cause 
elongation of the central isophotes would be either an error 
in frame alignment, or an unusual PSF from the AO cor- 
rection. However, both these effects can be dismissed: the 
globular cluster found I'.'lS from the galaxy centre is cir- 
cular and well fit by a 2D Gaussian of equal width in each 
dimension; no significant residual is seen after subtraction. 
We note that HST H-band images exists for NGC1399, but 
owing to a poorer spatial resolution (0'.'13) and PSF sam- 
pling, these structures are not observed. The non-spherical 
isophotes at the centre of NGC 1399 are a genuine feature of 
the galaxy. There appear to be peculiar kinematic features 
associated with the photometric anomaly: at r = —O'.'OS, the 
non-parametric VP contains substantial high velocity wings 
and at r = 0'.'14, the non-parametric VP is asymmetric with 
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an excess of receding velocity structure. We discuss possible 
explanations for these features below. 

Isop hote twists are c ommon in core-regions of core 
galaxies (iLauer et al.ll20()5ll . However, such twists are gen- 
erally smooth and seen on large scales, which is not the case 
here. The elongation could be a projection effect from the 
obscuration by a non-uniform distribution of dust. However, 
no obvious dust signature can be seen in either this Ks band 
image or archival HST V-band images. The density distri- 
bution of the dust would need to vary rapidly to cause such 
a sudden effect in the NIR on such small spatial scales (O'.'S 
or 29pc). However, neither of these explanations would give 
rise to the associated kinematic features. 

We have already seen that globular clusters (GCs) have 
been resolved near the centre of NGC 1399. Although the 
elongation of the nucleus does not appear to be separate 
to the central maximum, it is plausible that a GC happens 
to lie in the line of sight between ourselves and the galaxy 
centre and is unresolved from the central photometric max- 
imum of the galaxy. There would be a kinematic effect from 
such a chance alignment: one would expect a 'spike' in the 
VP at the systematic velocity of the GC caused by the low 
dispersion velocity structure of the globular cluster. How- 
ever, we see high velocity wings in the VP at r = — O'.'OS 
and a lopsided structure at r = 0'.'14 which is not the same. 

The elongation of the central isophotes is more likely to 
be a n offset centre as s een in a h andful of oth e r core galax- 
ies iLauer et al.ll20o3) . In fact, iLauer et alJ lf2005) argue 
that the offset centres detected in new WFPC2 p hotometry 
are a ll eccentric disks, analogous to that of M31 jTremaind 
^9^. Certainly, this would produce a strong kinematic sig- 
nature in the form of high velocity wings if our slit bisected 
the eccentric disk. Our slit does not completely bisect the 
peculiar photometry and nor does it coincide with the PA of 
the eccentric isophotes. However, the alignment of the AO 
reference star in the slit varied by around 1 pixel and this 
translates to a potential wandering of the slit, perpendicular 
to its length, by two image pixels. Thus, there is likely to be 
considerable contribution from the elongated photometry in 
the spectroscopic data. 

Keeping in mind the problems with registration between 
the image and kinematics, an eccentric disc at the centre of 
NGC 1399 would create high velocity wings at r = —O'.'OS, 
where the slit is closest to bisecting the anomaly. If the non- 
parameteric VP at this point is to be trusted, the velocity 
dispersion of just the wings (ignoring the central bulk) is 
in excess of 1000km s~^. Furthermore, if the stars in the 
eccentric disc rotate in the same direction, one would expect 
a lopsided VP at the pericentre of the ellipse if the disk was 
viewed near edge on. We do see a strong excess or 'hump' of 
structure receding away from the observer at r = 0'.'14 which 
is possibly the closest p ixel to the pericentr e of the elliptical 
photometry. Although, iLauer et al.l l|2005|) report no offset 
photometry in NGC 1399 with new WFPC2 photometry, 
the isophote ellipticity is seen to jump from approximately 
to 0.2 between 0'.'09 < r < O'.'l. 

We note that this eccentric disk hypothesis cannot ex- 
plain the decoupled kinematics discussed in Sec. 15.11 A co- 
herent disk can only survive well within the Sol of the BH, 
which we estimate to be ~ 0'.'3. 



6 DYNAMICAL MODELLING 

To obtain preliminary constraints on the mass of any BH in 
NGC 1399, we ignore the mild rotation gradient and elon- 
gated central isophotes and fit spherical dynamica l mod- 
els to our central kinematics combined with Graha m et alJ 
il998l) 's velocity dispersion profile. The latter extends to 70 
arcsec and provides important constraints on the galaxy's 
mass-to- light ratio. 

We assume that mass follows light, except at the galaxy 
centre where there can be a BH. The mass density distribu- 
tion is then 

p{r) = M.S{r)+rj{r). (16) 

Our goal is to find the range of BH masses M, and stel- 
lar mass-to-light ratios T that are consiste nt with our kine- 
matic s. We take j{r) from the models of iMaeorrian et alJ 
1^9^, which was obtained by deprojecting a composite sur- 
face brightness profile constructed by combining HST and 
ground-based photometry. Having this j{r) it is straight- 
forward to calculate the gravitational potential ipir) corre- 
sponding to 11611 for any choice of M, and T. Throughout 
this paper all mass-to-light ratios are in the V band. 

Our modelling procedure is a straightforward adapta- 
tion of the exte nde d Schwarzschild meth od described by 
iRix et alJ Jl997ft and lCretton et al.1 

(i) choose trial values for M. and T and calculate the 
corresponding potential tp{r); 

(ii) follow a representative sample of orbits in this ^p{r); 

(iii) find the weighted combination of orbits that min- 
imises the of the fit between the model and the observa- 
tions, subject to the constraint that each orbit carries non- 
negative weight; 

(iv) assign the likelihood exp(— ^x^) to the potential ip- 

We do not impose any regularization in the third step. By 
considering a range of plausible trial potentials and compar- 
ing their relative likelihoods, we obtain constraints on both 
M. and T. 

6.1 Orbit distribution 

Apart from the parameters (M. , T) describing the potential, 
the other unknown in our models is the distribution func- 
tion (DF) /(x,v), defined such that /(x, v)d^xd^v is the 
luminosity of stars in some small volume d^xd^v of phase 
space. Jeans' theorem tells us that if our models are to be in 
equilibrium their DFs can depend on (x, v) only through the 
integrals of motion £ and J, the binding energy and angular 
momentum per unit mass. We make the stronger assump- 
tion that the local stellar velocity distribution is symmetric 
about the Vr axis, so that the DF depends on J only through 
its magnitude. We then discretize the DF as a sum of delta 
functions, 

f{£, j') = J2J2 ^^^^^^ ~ - "^So, (17) 

on a regular ne x nj grid in phase space. The points £i are 
chosen through £i = tp{ri) with the spaced logarithmically 
between 1 pc and 100 kpc. This ensures that orbit apocen- 
tres are approximately uniformly distributed among each 
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decade in radius. For each £i there are nj values of angular 
momentum, with J^j running linearly between and J^{£i), 
the angular momentum of a circular orbit of energy £i. To 
avoid a rash of indices we also write the double sum 1171 as 
a single sum over n = ne x nj points: 



(18) 



6.2 Observables 

Having a trial potential tpir) and a set of DF components 
{£k, Jk)y we calculate the unnormalised, psf-convolved VP 
histogram of each component at the projected radius Ri of 
each of our riAO ~ 31 kinematical data points: 



r(fe) 



I' 



dvz / dxdypsf{Ri — X, —y) 



X / dVxdvyS{£ — £k)S{J^ — Ji). 



(19) 



Here we use a rectangular co-ordinate system {x,y,z) with 
origin O at the galaxy centre and Oz-axis parallel to lines 
of sight. For our standard models each histogram has n„ = 



24 bins of width Vj 



Vj-l 



50km s ^, with innermost 



bin edge at vi = 0km s~^, the galaxy's systemic velocity, 
and outermost edge at 1125 — 1200km s~^. The function 
psf(Aa;, Ay) is the two-dimensional off-source point-spread 
function of Sec. 12.51 and Table Armed with the L^j , the 
unnormalised VP histogram of a model with DF (/i, . . . , /„) 
is simply 



(20) 



The normalisation constant is the psf-convolved surface 
brightness. 



/(i?o = /. = E^^^*'' 



(21) 



6.3 Fitting models to observations 

6.3.1 Gauss-Hermite coefficients 

Calculating the Gauss-Hermite coefficients of our models is 
simple: the (unnormalised) contribution of the fc"^ DF com- 
ponent to the i*'^ VP is given by i'^^ with j = 1, . . . , n„, 
from which equation IIH allows us to calculate this com- 
ponent's contribution /i'*^' to {hj}. However, as we have 
explained in section 13.11 the observed Gauss-Hermite co- 
efficients are not independent. It takes only a little effort to 
include the effects of the covariances among the {hj} into 
the modelling. 

Let us define the column vector h = (ho, . . . , hN)'^ and 
let h be the vector of coefficients that minimise the Xs of 
equation ((7||. For fixed normalisation k and continuum pa- 
rameters Ci, this Xs is a quadratic form in the hi: 



x'(h) 



2 

Xmin 



AT 

+ ^ VA4e,-(h-h)]2 
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i=0 



(23) 



d^Xa /9hidhj has eigenval- 



where the Hessian Mij 
ues Ao, . . . , Ajv with corresponding eigenvectors eo, . . . ejv. 
Therefore, the new parameters 

h'^ =e, h (24) 



have independent errors Aj = y^2/Xj. 

Summing the results from I23II for each of our measured 
VPs, 



Xgh = 



"AO -'V 

EE 



A. 



(25) 



where we have dropped the Xmin terms and have introduced 
the h'^^\ which are related to the ft^^' by (I24I I. The model's 
predicted VPs in 12511 are normalised by the "observed" local 
surface brightnesses U. The latter are difficult to extract 
from spectroscopic observations. Instead, we obtain them by 
convolving the photometric profile used to obtain our j(r) 
profile with the effective spectroscopic PSF. Of course, a 
reasonable model must also fit this I{R), so we add 



where, by analogy with 119II . 



dxdypsi{R,~x,-^y) / d\S{£~£k)S{j''-Jl).(22) 



r(fc) 



We explain how we evaluate the multiple integrals 1191 
and I22II in iMagorrian et aH i2005l) . The spherical symme- 
try of our models means that we can afford to calculate 
these projection coefficients for each of our 31 observed 
radii directly. More sophisticated axisymme tric models (e.g., 
JCretton et alJll999l: FGebhardt et al.ll2003l) 'l usually resort 
to introducing subgrids to store intermediate quantities in 
this calculation. 

Our treatment of Graham et al's velocity dispersion pro- 
file is similar. We assume that each of their data points 
measures the second moment of the VP convolved with a 
Gaussian PSF with FWHM 2 arcsec. The (unnormalised) 
second moment of each DF component is given by 12211 with 
an extra factor of v1 inside the innermost integral. 



2 
Xi 




(fe) 



A/. 



(26) 



to 1251 1. somewhat arbitrarily assigning errors A/i — 10"'' /i. 
Finally, we add one more term, xgi to measure how well our 
models fit both the unnormalised second moments from Gra- 
ham et al.'s data and the corresponding psf-convolved sur- 
face brightnesses, again assuming fractional errors of 10~^ 
in the latter. Adding all these measurements together, the 
final x^ of our model. 



Xm[^,/fc] = xio + Xl +XG, 



(27) 



which, for fixed potential 1/), depends quadratically on the 
orbit weights /fc. We use a standard non-negative least- 
squares algorithm iLawson fc Hanson 1974') to find the non- 
negative set of /fe that minimise it. Unlike most other orbit- 
superposition methods, we do not include the luminosity 
density j(r) in this fit explicitly, but use the surface bright- 
ness profile I{R) instead. Our reason for this is that, un- 
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like j{r), I{R) is (in principle) directly measurable and can 
therefore be assigned meaningful error bars. We do not ex- 
pect real galaxies to have perfectly constant mass-to-light 
ratios and, in the absence of anything better, use j{r) merely 
to estimate the stellar contribution to the overall potential. 



6.3.2 VP histograms 

Unlike the Gauss-Hermite coefficients above, there is no 
simple way to compare the VP histograms found in 
against models: in addition to the unavoidable correlations 
among the Li, there are complicated biases introduced by 
the penalty function HSU . So, we simply refit the VP his- 
tograms using the same riv velocity bins for which we cal- 
culate the models' VPs in H6.2I reflecting the VPs about 
V — in order to make them symmetric. The procedure is 
as follows: 

(i) First find the best-fit template fraction, continuum 
level and normalisation: find the /, d and smoothed Li that 
minimise the penalised Xp = Xs + -P[^i] for A; = 1, where 
Xs and the penalty function P are given by equations Q 
and list respectively. The best-fit normalisation factor is 
then = ~ Vi)Lt. 

(ii) Holding k, the d and / fixed, find the formal best-fit 
, Ln^) to the unpenalised Xs 0- 



VP histogram (Li 



The best-fit VP L is in general unphysical with many Li < 
0, but we use it only because it locates the minimum of 
the Xs quadratic form. Let us write our formal best-fit VP 
histogram as the column vector L = (Li )"^ and 
consider another L = (Li, . . . ,L„„)"^. Once fe, / and the d 
have been fixed, equation O becomes 



Xs(L) = Xmi 



+ 



L — Gi ■ L 



(28) 



where are the eigenvectors of the Hessian d^x^ / dLidLj 
and the eigenerrors Ai are related to the eigenvalues Xi 
through Ai = yj2/Xi. So by taking the projections 



L'. 



(29) 



of L along the full set of "eigenVPs" e^, one test directly 
how well it reproduces the observed galaxy spectrum. We 
note that our method is essentially a restatement of the 
work of^Do Riicko & Doionghc ( 1998) using the language of 
eigenVPs. Figure|S|plots the first few eigenVPs of one of our 
spectra. Unlike the terms in a Gauss-Hermite expansion @, 
they do not taper off rapidly at high velocities. 

Summing the results of H28|l for each VP and neglecting 
the Xmin terms, a model with DF (/i, . . . , /„) has 
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(30) 



where iij*^' is obtained from L^*^' through 12911 . Apart from 
replacing Xgh in equation 1271 by 13011 . our procedure for fit- 
ting models to VP histograms is identical to that for Gauss- 
Hermite coefficients. 




Figure 6. Plot of the first four eigenVPs of one of our spectra 
obtained by diagonaUzing Xb of eq. J7J. Any (symmetric) VP can 
be expressed as a weighted sum of eigenVPs feg. 1^91 . in which 
case the errors in the weights are independent. We have divided 
each of the VPs plotted here by its corresponding eigenerror so 
that the scale of each gives a direct indication of the uncertainty 
in its weight. 



6.4 Results 



We have calculated the projection coefficients L['^\jp) and 

ll''\ip) for ne xnj = 200 x 40 DF components (eq.ITTJ in a 
range of potentials tp. Our main results below are obtained 
using coefficients calculated for BH masses M./10^Mq = 
0, 0.25, 3.75 with a single mass-to-light ratio. To — 
8.5T0. Since the projection coefficients scale straightfor- 
wardly with mass, we can use the method of Hti.3l to fit mod- 
els with other values of T provided we remember to scale 
Af, by T/To and the bins vj of the velocity histograms 1191 
by y/rjTo. 

On figure [7| we plot the result of using 1271 to fit 
the Gauss-Hermite coefficients {ho, h2, . . . , /ijv) for the cases 
A*' = 4 and A = 6. In each case we have extracted VPs 
from spectra assuming that hi = for i > N. Marginal- 
izing T, the former yields a moderately firm BH mass 
M. ~ 2.3to l X 10® Mq. Including he in the fit, however, 
results only in the upper bound. A/. < 1.6 x l(f Mq, sta- 
tistically consistent with the A = 4 result. To test whether 
this could be a result of the best-fit template fraction /, 
continuum d or normalisation k (Sec. ^ changing as we 
change A, we have repeated our A = 6 model fit omitting 
the coefficient h'j (eq. 1241 with the largest eigenerror Aj . 
The resulting Xm(A/. ,T) contours are broadly similar to 
the A = 4 case, suggesting that the shift in the location of 
the minimum in Xm is a genuine effect caused by the extra 
information contained in the fee parameter, and not by sys- 
tematic changes in the fitted continuum or normalisation. 

These results are perhaps unsurprising when one re- 
calls that our results in Sec. 2] indicate that NGC 1399 
has very strongly non-Gaussian central VPs. Although it 
would be a mildly interesting exercise to test the effects of 
fitting to hs and even higher-order terms, we now simply 
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Figure 7. XmC^*''"") contours (eg. I27i obtained by fitt ing dynamical model s to Gauss— Hermite parametrisations fSec. I^!T1 of our 
spectra and to the velocity dispersion profile measured bv lGraham et all ^^9^. Successive contour levels have = 1- The models 

have 200 X 40 DF components I17i and fit {ho, h2, ■ ■ ■ , /ijv) extracted from spectra under the assumption that hi = for i > N. The left 
panel shows results for N = A, the right for N = 6. The x-axis is both cases is the BH mass scaled by Tg.s = T/8.5Tq. Mass-to-light 
ratios T are V band. 
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Figure 9. Anisotropy parameter /3 = l — Ug/a^ of a typical model 
(M, = 10^ Mq, T = 8.5T0), averaged over five shells per decade 
in radius. 

drop the Gauss-Hermite parametrisation and turn to fit- 
ting eigenVPs. Fitting models to the full set of eigenVPs 
yields the result plotted in the bottom-right of fig. |H| The 
best-fit BH mass, 1.2^Q g x IO'^ATq, is consistent with the 
results from the Gauss-Hermite fits, but has smaller error 
bars. The other panels on the figure show the effect of rear- 
ranging the eigenVPs in order of increasing eigenerror and 
fitting only the first for each spectrum. The case = 3 
yields results that look qualitatively similar to our sixth- 
order Gauss-Hermite fit, while adding one more eigenVP 
introduces a useful lower-bound on the BH mass. The re- 
sults for N = 5 are very similar to fitting the full A'^ = 24: 
in fact, it is impossible to distinguish between = 6 and 
= 24 by eye. 

Since the the A'^ = 24 model is essentially a direct 
fit to the galaxy spectrum we adopt its best-fit M, = 
l-2lo g ^ our best estimate of the BH mass in 

NGC 1399. Figure 1^ plots the anisotropy parameter of one 
of the best-fit models in this range. Our models have mod- 



erate radial anisotropy {f3 « 0.3) be tween 2 arcsec and 30 
arcsec, similar to the results found bv lSaelia et al.l i200Cfl in 
their models of NGC 1399. At larger radii the orbit distri- 
bution in our models becomes tangentially biased, which is 
what one expects from fitting a constant mass-to-light ra- 
tio model to a galaxy with a massive dark halo. Much more 
surprisingly, however, our models become extremely tangen- 
tially biased in the innermost arcsec, which is where we find 
the interesting kinematic and photometric features. 



6.5 Tests 

Our best-fit BH mass remains unchanged if we change the 
number of DF components we use in the models: models us- 
ing 100 X 20 components instead of 200 x 40 yield the same 
result. We find no evidence for the flat-bottomed pro- 
files claimed bv lValluri. Merritt. fc EmsellemI ||2004|) (albeit 
for the axisymmetric case), even though our models sample 
{£, J^) phase space a factor ~ 20 more densely than theirs. 
We have tried varying the extent and widths of the velocity 
bins we use when we fit VPs, but the BH mass is unchanged 
whether we take 24 fine 25km s~^-wide bins extending to 
1200km s"^or 20 coarse 100km s ^-wide bins extending to 
2000km s"\ 

The assumed width of the velocity bins does affect the 
models in a slightly more subtle way, however. When we 
rescale our base To — 8.5T0 models to a new mass-t o- light 
ratio T, we scale their velocity bins by an amount ^JT/Tq 
and re-extract VPs from spectra using the the new bin 
widths. Therefore, values of Xs(L) from 12811 for different T 
are not strictly directly comparable. To test the effects of 
this we have calculated coarser 100 x 20-component models 
using 24 bins of fixed width 50km s^^on a grid of poten- 
tials with T/Tq running from 7.2 to 10.6 in 0.2 steps and 
M./10^Mq running from to 3.75 in 0.25 steps. The only 
notable difference between the resulting XmiM,, T) and the 
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Figure 8. As for fig. El but fitting to eigenVPs (' il6.:i2l instead of Gauss-Hcrmitc coefficients. Tfie first tfiree panels show the results of 
fitting only to the first 3, 4 and 5 eigenVPs of each spectrum. The last shows the results of fitting to the full spectrum by using all 24 
eigenVPs. 



plots in fig.|H]is that the range of acceptable T increases to 
(9 ± 1)T0. The BH mass is unaflFected. 



6.6 Caveats on BH mass 

Despite these reassurances, there nevertheless are some 
shortcomings of the models and data upon which this BH 
mass is based: 

(i) We follow the usual "extended-Schwarzschild" proce- 
dure from which the majority of existing stellar-dynamical 
BH masses have been obtained and consider only the very 
best orbit distribution fk for each potential. This best-fit 
distribution is typically very spiky, with only ~ 140 non- 
zero fk out of 200 X 40! One way around th is would be to 
apply some kind of regularization to the fk iThomas et alJ 
I2OO3), but the b iases introduced by this procedure are not 
weU understood (IVaUuri. Merritt. fc Emsellemlliooi) . 

(ii) Because of the huge freedom in fitting the fk, the 
best-fit model should have a very low Xm- For example. 



in Monte Carlo experiments with synthetic datasets of toy 
galaxy models we typically find xL ~ A'^data/3, where A''data 
is the number of data points fit in the models. For our real 
NGC 1399 data, however, our best-fit xgh s^nd Xh usu- 
ally at least as big as the number of parameters we use to 
describe our kinematics. This is probably due to our neglect 
of the systematic errors (Sec. l2.3l and l4.2ll . In contrast, the fit 
to the full surface brightness profile and to the outer disper- 
sion profile is astonishingly good, with xI + Xg ~ 3, making 
our total Xm relatively low. 

(iii) Finally, our models assume that galaxy is spherical 
and non-rotating, despite the clear evidence to the contrary. 
Nevertheless, the models do seem to require a strong bias 
towards circular orbits in the central 0.5 arcsec. 



7 CONCLUSIONS 

Using NAOS-CONICA at the VLT, we have successfully 
measured the central kinematics within the Sol of NGC 1399 
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(r ~ 0'.'34) with a resolution (FWHM) of ~0'.'15 (14pc, Sec 
12.51 using adaptive optics correction on a bright reference 
star l/.'S away. 

Alone, the kinematics extracted from the Cal feature at 
and the CO bands after 2.3/im establish the pres- 
ence of velocity gradient within a radius of ~0'.'5 (48pc), 
suggestive of a kinematically decoupled core. 

Ks band imaging reveals offset and elongated isophotes 
within a radius of 0'.'2 (19pc) that are not visible in H- 
band HST images. Such inner structure is reminiscent of 
that seen in other core ellipticals. The non-parametric VPs 
corresponding to this region also show an unusual velocity 
structure that may be consistent with the presence of an 
eccentric disk around the BH, akin to that of M31. 

We have demonstrated that errors in the Gauss- 
Hermite coefficients hj extracted from real galaxy spectra 
are not independent, and have shown a simple way of taking 
the covariances among the hj into account when fitting dy- 
namical models. The VPs near the nucleus of NGC 1399 are 
strongly non-Gaussian, however, and are not well described 
by a low-order Gauss-Hermite expansion. We show that the 
"eigenVPs" obtained by diagonalizing O are a more useful 
way of describing VPs, at least for numerical purposes. 

Subject to the caveats of §6.6, our best estimate for 
the mass of the BH in NGC 1399 is 1.2^°;^ x 10^ Mq, ob- 
tained by fitting spherical dynamical models directly to our 
observed spectra. The models are based on the usual exten- 
sion of Schwarzschild's method to the problem of potential 
estimation. Taken at face value, they place this galaxy on 
the M.-o" plane mid- way between the predictions of T02 and 
FF05, being consistent with both. The best-fit model also 
becomes extremely tangentially anisotropic in the innermost 
C.'S. 

We have demonstrated that AO observations are a 
viable alternative to IfST to when measuring black-hole 
masses and can break the mass-anisotropy degeneracy even 
in the most massive, non-rotating elliptical galaxies. The 
difficulty in interpreting these long-slit data emphasises the 
need for high SNR AO assisted integral-field observations to 
further understand the kinematic and photometric features 
discovered here. 
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